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ABSTRACT 

We report here the first systematic study of the supercluster-void network in the 
ACDM concordance cosmology in which voids and superclusters are treated on an 
equal footing. Superclusters are defined as individual members of an over-dense excur- 
sion set and voids are defined as individual members of a complementary under-dense 
excursion set at the same density threshold. We determine the geometric, topological 
and morphological properties of the cosmic web at a large set of density levels by 
computing Minkow ski functionals for every supercluster and void using SURFGEN 
ijSheth et al.ll200^ . The properties of the largest (percolating) supercluster and the 
complementary void are found to be very different from properties of individual su- 
perclusters and voids. Individual superclusters totally occupy no more than about 5% 
of the total volume and contain no more than 20% of mass if the largest supercluster 
is excluded. Likewise, individual voids totally occupy no more than 14% of volume 
and contain no more than 4% of mass if the largest void is excluded. Although super- 
clusters are more massive and voids are more voluminous the difference in maximum 
volumes is not greater than by an order of magnitude. The genus value of individual 
superclusters can be ~ 5 while the genus of individual voids can reach ~ 40, implying 
significant amount of substructure in superclusters and especially in voids. One of our 
main results is that large voids, as defined through the density field (read dark matter 
distribution) can be distinctly non-spherical. 
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1 INTRODUCTION 

One of the great observational discoveries of recent times is 
the realization that we live on a cosmic web which is embed- 
ded in an accelerating Universe. Although an accelerating 
universe was o riginally established using high red shift type 
la supernovae fenop et al.l2003tlTonrv et alj20o3) . the case 
for a 'dark energy' dominated universe has received inde- 
pendent sup port from observatio ns of the cosmic microwave 
background iSoergel et aLllioOSf) combined with studies of 
galaxy clust ering in the two degre e field (2dF) galaxy red- 
shift survey iEfstathiou et al1l2002h . Although the nature of 
dark energy is still an open issue, a cosmological constant ap- 
pears to agree very well with all current observations which 
indicate f^A — .73, Q.m — 0.23 and Q.b — 0.0 4 (for reviews of 
dark energy see lSahni fc Starobinsk^1 ll2000D : ISahnl J2002D 1. 
Fully three dimensional large scale galaxy catalogues reveal 
that the cosmic web consists of an interpenetrating network 



of superclusters and voids. It therefore becomes important 
to understand and quantify the geometrical and topological 
properties of large scale structure in a ACDM cosmology in 
a deep and integrated manner. 

The main aim of the present study is to study the 
supercluster-void network in ACDM cosmology with empha- 
sis on the sizes, shapes and topologies of individual super- 
clusters and voids. We shall also study the percolation prop- 
erties of the full excursion set sampling the entire density 
field and quantify our results in terms of Minkowski func- 
tionals (hereafter, MFs). A study such as the present one 
will help us in comparing theory against observations. It 
will also help us to distinguish the salient features of ACDM 
cosmology from sister cosmologies in which dark energy is a 
function of time. (These issues will be discussed in detail in 
companion papers.) 

It may be appropriate to mention that this is the first 
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comprehensive analysis of large scale structure geometry and 
morphology in which over-dense (superclusters) and under- 
dense (voids) regions are treated on a completely equal foot- 
ing. Earlier studies have emphasized either over-densities 
(clusters, superclusters) or under-densities (voids) as a result 
of which the methods used for the analysis of these two com- 
plementary entities (superclusters/voids) often vary greatly 
in the literature. Thus over-dense regions have been studied 
using correlation functions, minimal spanning trees, shape 
functions etc., whereas voids have been constructed from 
point processes using elaborate boundary and volu me filling 
techni ques (see e.g. Sahni & Coles (1995); Martinez fc Saarl 
i200lf) and references therein) . Although the above methods 
do provide us with some insights into supercluster-void mor- 
phology, we shall follow an alternative route in this paper 
and study the properties of the supercluster-void network 
using the density field as our main starting point and fun- 
damental physical quantity. The reason for this is two fold, 
firstly, density fields can be easily constructed from point 
data sets using, for example, cloud in cell (CIC) techniques 
(this is true both for data from N-body simulations, which 
we examine in this paper, as well as for galaxy distributions 
in three dimensions). Secondly, an elaborate surface mod- 
eling scheme, SURFGEN (short for 'surface generator'), al- 
lows us to determine the geometry, morphology and topol- 
ogy of excursion sets defi ned on a density field in a very 
comprehensive manner fiSheth et aLlEoOSl . Applying SUR- 
FGEN to the density distribution in the ACDM model al- 
lows us to develop deep insights into the distribution of large 
scale structure in this model, which can be quantified using 
MFs and Shapefinders. Working with the density field also 
permits us to determine the morphological properties of the 
full excursion set describing the supercluster-void network. 
More detailed information is then gleaned at one particu- 
lar threshold (usually associated with percolation) at which 
shapes and sizes of individual superclusters and voids yield 
rich information about properties of the cosmic web to which 
we belong. 

Concretely, we study the large scale structure of the uni- 
verse by considering the geometry and topology of isodensity 
surfaces 5(x) = 5p(x)/p = const. At a given threshold 5th 
regions having higher than threshold density (5 > 5th) will 
be called "superclusters", while regions with 5 < <5th will 
be called "voids" . Thus, we define superclusters and voids as 
over-dense and under-dense connected regions bounded by 
one (or several) surfaces of constant density. This definition 
broadly corresponds to other definitions of superclusters and 
voids used in the literature but differs in details. Apart from 
obvious differences with superclusters and voids of galaxies 
in redshift space our approach specifies neither a particular 
density threshold nor the shapes of the structures. Despite 
these differences we call the objects of our study superclus- 
ters and voids mostly because they are nonlinear structures 
having sizes, volumes and masses roughly corresponding to 
superclusters and voids of galaxies. 

We study superclusters and voids at a large number of 
density thresholds and construct the isodensity surfaces at 
every thresh old to best accuracy. In c o ntrast to many stud- 
ies (se e e.g. iBlumenthal et all lll992h : IColdberg fc Voeelevl 
l|2003l) l. we do not 'cook up' voids or superclusters with 
predefined shapes but isolate individual objects from the 
dark matter density field obtained in the N-body simula- 



tions by building the isodensity surfaces. Apart from uni- 
formly smoothing the density field with a Gaussian filter, we 
do not introduce any factors that may affect the shapes or 
substructure in superclusters and voids of the cosmic density 
field. Filtering the density field may erase some small-scale 
features but it certainly does not introduce any new struc- 
tures. Thus, we know beforehand that the real structure 
can be only richer and more elaborate compared to what 
we study after smoothing. Filtering high frequency modes is 
virtually implied in every physics study. This approach can 
be also viewed as an application of the standard excursion 
set technique to non-Gaussian three-dimensional fields. 

One of the goals of this work is a quantitative and 
objective testing of some stereotypes routinely used in 
cosmology. Here are a few examples: both N-body simula- 
tions and redshift surveys are characterized by filamentary 
and sheet-like or pancake-like structures; voids are quasi- 
spherical bubbles in the density and/or galaxy distributions; 
and expand faster than the universe as whole; voids occupy 
most of the volume in the universe. The first problem 
arises when one tries to address these issues in the absence 
of conventional definitions. Frequently the context of the 
origin and usage of a particular assertion is either obscure 
or forgotten. Proving that some of the above cliches cannot 
be true does not require elaborate N-body simulations 
or sophisticated analysis. For instance, voids cannot both 
occupy most of the universe and expand faster than the 
universe because this would require the expansion of voids 
in comoving space. But this is impossible since the comoving 
volume is conserved ! Therefore, the expansion of some voids 
must clearly be at the expense of some of the others. This is 
one of the trivial (but imp ortant) conclusions arising from 
the adhesion model l|&irbato\^^ai chev fc S handariiil 

Il98^ |Shandarin & ZeldovichI 119891 : 

Sah ni. Sathvaprakash & Shandarin 199^. At the same 
time it is not necessary for the interior of the squeezed 
voids to shrink. The interior of the void may continue 
to expand faster than the Hubble fiow but the void 
eventually vanishes because its boundaries move inwards 
I Sahni. Sat hvaprakash & Shandarin 1994 . ■1995 a b). So 
as the void collapses it is likely to become increasingly 
non-spherical. Thus, a single expanding void model while 
being reasonable for a crude estimate of substructure 
growth may be completely misleading for estimating the 
sizes of a realistic ensemble of voids. 

Pancakes or sheet-like structures remain an unsolved 
controversy since their theoretical prediction by Zel'dovich 
in 1970. On the one hand there were numerous asser- 
tions that they had been detected in the N-body simu- 
lations or redshift surveys most of which were based on 
visual impressions. On the other, many believe that pan- 
cakes are not clearly seen in any realistic N-body simula- 
tions and have not been detected objectively. It is worth 
noting that the simulation of the structure in the hot 
dark matter model did not show the existence of pancakes 
jKlvpin fc Shandarinll98il which was interpreted as the re- 
sult of weak singularities corresponding to pancakes com- 
pared to con siderably stronger singularities correspo nding 
to filaments llArnol'd. Shandarin fc Zel'dovichl Il982h . The 
coarse graining unavoidable in N-body simulations erases 
pancakes while filament s survive. This argument has been 
essentially repeated by iBond. Kofman fc PoeosvarJ (Il996l) 
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and elaborated for the probabilistic interpretation of the 
large-scale structure. In this study of the ACDM model we 
do not see pancakes either. 

The large-scale structure in the dark matter density 
field is certainly different from the one observed in galaxy 
redshift surveys. Galaxies are at best pointwise tracers of 
the parent continuous dark matter density field. It is likely 
that they are biased tracers, and that they do not display 
real physical structures but structures which are strongly 
distorted by the mapping into redshift space. Thus, the ob- 
served superclusters and voids are not physical but only 
apparent objects as the retarded motion of some planets 
is only apparent but not real motion. At present there are 
only two methods of investigating real large-scale objects: 
by reconstructing the real space density field from peculiar 
velocities of galaxies and by investigating the density fields 
in the models that are consistent with observations such as 
ACDM. 

We consider this study of the dark matter density field 
in real space as a first step in a systematic study of the 
morphology of large-scale structure. It will be followed by 
studies of mock and real galaxy catalogues. Understanding 
the morphology of the dark matter distribution in real space 
is an important component in understanding the physical 
processes determining the formation of galaxies and their 
motions as well as for building theoretical models of super- 
clusters and voids. 

Due to the fairly large smoothing scale adopted here, 
the over-density in superclusters ranges from 5 ~ 1 to 5 ~ 
10, which makes them more extended ( few Mpc) and 
considerably less dense than galaxy clusters in (for instance) 
the Abell catalogue. 

The rest of the paper is organized as follows. In Section 
|5| we briefly describe the density fields used in this study. 
In Section |3 we define and very briefly discuss the mor- 
phological measures which we apply to study the cosmic 
web. Sec.|l]provides the relation of the filling factors used in 
this study with the probability density function. Sec. 5 dis- 
cusses the morphological properties of individual superclus- 
ters and voids. Sec. 6 describes substructure in superclusters 
and voids. A discussion of some properties of Minkowski 
functionals is contained in Sec. 7. Section 8 summarizes our 
main results. 



2 DARK MATTER DISTRIBUTION IN VIRGO 
SIMULATIONS 

We use dark matter distributions in a flat model with flo — 
0.3, Qa = 0.7, h = 0.7 (ACDM ). The initial spectrum was 
taken in the form suggested bv lEfstathiou. Bond fc White! 
il992h with the shape parameter F = 0.21. 

The amplitude {as = 0.9) of the power spectrum in the 
model is set so as to reproduce the observed abundance of 
rich galaxy clusters at the present epoch. A detailed discus- 
sion of the cosmological parameters and simulations can be 
found in Jenkins ct al. (199^. 

The density fields generated by the VIR GO simulations 
have been studied in several papers (see e.g . ISpringel et alJ 
Jl998h :ls coccimarro. Couchman fc FriemarJ (119991) 1: this is 
the first systematic study of global and local morphology 



when both over-dense and under-dense regions are measured 
by the same method. 

SURFGEN operates on three-dimensional pixelized 
maps. Therefore we first generate the density field from the 
distribution of dark matter part i cles. T his process was de- 
scribed in detail in ISheth et alj JioO^); here we present a 
brief summary. The data consist of 256^ particles in a box 
of size 239.5 fe"^Mpc. We fit a 128^ grid to the box. Thus, 
the size of each cell is 1.875 /t~^Mpc. Here we fol low the 
smoothing technique used by ISpringel et al.l lll99^ which 
they adopted for their preliminary topological analysis of 
the Virgo simulations. In the first, we apply a Cloud in Cell 
(CIC) technique to construct a density field on the grid. 
Next we smooth this field with a Gaussian kernel which of- 
fers us an extra smoothing length-scale. We study the field 
smoothed with Ls = 5 /i~^Mpc but also present the global 
MFs for the field smoothed with 10 /i~^Mpc for compar- 
ison. The scale of Ls = 5 /i~^Mpc is a fiducial smooth- 
ing scale in many studies of both density fields in N-body 
simulatio ns and galaxy fie l ds fro m redshift surveys; see for 
example, iGrogin fc Gelled ||2000D . The Gaussian kernel for 
smoothing that we adopt here is 

Since the kernel is isotropic, it is likely to diminish the 
true extent of anisotropy in filaments and pancakes. This 
effect could be minimized by considering anisotropic kernels 
and/or smoothing techniques based on the wavelet trans- 
form. An even more ambitious approach is to reconstruct 
density field s using Delaunay t essela tions using a technique 
reported by Ivan de Wevgaer'3 tOOj) . Density fields recon- 
structed in this manner appear to preserve anisotropic fea- 
tures and may therefore have some advantage over con- 
ventional 'c loud-in-ceir techniques fol lowed by an isotropic 
smoothing JSchaap fc WevgaerdEoOoll . Thus, as far as one 
of the goals is to utilize the geometry of the patterns to 
discriminate between the models, such smoothing schemes 
should prove more powerful. We hope to return to these 
issues in further publications. 

We scan the density fields at 99 values of the density 
threshold, all equispaced in the filling factor FFc defined as 

FFc{Stii) = -^ I e{S^5Tu)d\, (2) 

Vtot J 

where Q{x) is the Heaviside Theta function and 5 = {p — 
p)/p is the density contrast. The supercluster filling factor 
FFc measures the volume-fraction in regions which satisfy 
the 'supercluster' criterion S^c ^ 5th. In the following, we 
use FFc along with the void filling factor FFy = 1 — FFc 
as a parameter to label the density contours. The relation 
between FFc and density contrast threshold can be seen in 
Fig.IH 

At each level of the density field (labeled by FFc or 
FFv), we construct a catalogue of superclusters (over-dense 
regions) and voids (under-dense regions) based on a grid 
realization of the Friends of Friends (FOF) algorithm. 

Next we (i) run the SURFGEN code on each of these 
superclusters/ voids to model surfaces for each of them and 
(ii) determine the Minkowski Functionals (MFs) for all su- 
perclusters/voids at the given threshold (these are referred 
to in the literature as partial MFs). Global MFs are partial 
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MFs summed over all superclusters. Thus, at each level of 
the density, we first compute the partial MFs and then the 
global MFs. 



3 MORPHOLOGICAL PARAMETERS 

Since the 1970s, theory, N-body simulations, and, most im- 
portantly, galaxy redshift surveys have strongly suggested 
that the components of the large-scale structure can be 
roughly placed in three classes: compact quasi-spherical or 
slightly elliptical structures like Abell clusters, long fila- 
ments like the famous 'bridge' connecting the Coma cluster 
and A1367, (Gregory fc Thoinoson (1978)) and voids. There 
have also been claims that pancake-like concentrations of 
galaxi es have been observed iFairalJ |l99i);|MartiHizASii3 
(|2001j). The voids have often been claimed to have quasi- 
spherical or slightly elliptical shapes. Most of these claims 
have been based on visual impressions. In particular, statis- 
tics used in studies of voids often assumed that voids are 
spherical or close to spherical thereby precluding any other 
possibility. 

In this work we study the geometry and topology of 
the regions bounded by the isodensity surfaces and there- 
fore make no prior assumptions about the shapes of super- 
clusters and voids. It is worth noting that some regions may 
have more than one boundary surface and possess nontrivial 
topology of the boundaries. The complete characterization 
of an arbitrarily complex region in three dimensions obvi- 
ously cannot be achieved if only a few numbers are used. At 
best one can try to design some basic characteristics that 
serve a particular purpose. Our purpose is to provide basic 
measures suitable for quantification of typical components 
of the large-scale structure; superclusters and voids. 

Four Minkowski functionals are effective non- 
parametric descriptors of the morphological properties 
of surfaces in three dimensions iM cckc. Buchcrt & Wagner. 
ll994l:lMatsubarall2003l : ISheth et al]l200A . They are 

• Volume V enclosed by the surface, S, 

• Area A of the surface, 

• Integrated mean curvature C of the surface, 



da. 



(3) 



where Ri and R2 are the principal radii of curvature at a 
given point on the surface. 
• the Euler characteristic 



1 

2^ 



{R1R2) 



da. 



(4) 



Although for our purpose the Euler characteristic is a more 
convenient quantity than genus, G (see Appendix for a dis- 
cussion) the latter has been used more often in cosmology. 
The genus is uniquely related to the Euler characteristic 
G = 1 — x/2 and thus carries exactly the same informa- 
tion. We measure the above parameters for every region in 
both over-dense and under-dense excursion sets at 99 density 
thresholds equispaced in the filling factor from FF — 0.01 
to FF = 0.99. 

A s demonstrated in |SalinL_Sathvaprakash,_£ ^handarinl 
il998h : ISathvaprakash. Sajm^^ handarinT i^g^ particu- 
lar ratios of Minkowski functionals called "Shapefinders" 



provide us with a set of non-parametric measures of sizes 
and shapes of objects. Therefore, in addition to determining 
MFs we shall also derive the Shapefinders, T (Thickness), 
B (Breadth) and L (Length) defined as follows: 

The three Shapefinders describing an individual region 
bounded by one or several isolated surfaces of constant den- 
sity have dimensions of length and provide us with an esti- 
mate of the regions 'extension': T is the shortest and thus 
describes the characteristic thickness of the region or ob- 
ject, L is typically the longest and characterizes the length 
of the object; B is intermediate and can be associated with 
the breadth of the object. This simple interpretation is ob- 
viously relevant only for fairly simple shapes. The choice of 
the coefficients in eq. |^ results in a sphere having all three 
sizes equal to its radius T = B = L = R. A triaxial ellipsoid 
has values of T, B and L close but not equal to the lengths 
of its three principal semi-axes: shortest, intermediate and 
the longest respectively. It is worth noting that T, B and L 
are only the estimates of three basic sizes (semi-axes) of an 
object which wor k quite well on such objects as a triaxial el- 
lipsoid and torus l|SaImij_SatlixaErak^^ 
ISathvaprakash. Sahni fc Shandarinlll998HSheth et al.ll200a) 
but no three numbers can describe an arbitrary, complex 
three-dimensional shape. 

An indicator of 'shape' is provided by a pair of dimen- 
sionless Shapefinder statistics 



B-T 
B + T' 



F : 



B 



L + B' 



(6) 



where P and F are measures of Planarity and Filamen- 
tarity respectively {P,F 5C 1). A sphere has P = _F = 0, 
an ideal filament has P = Q, F = 1 while P = 1,_F = 
for an ideal pancake. Other interesting shapes include 'rib- 
bons' for which P F ^ 1. When combined with the genus 
measure, the triplet {P, F, G} provides an example of shape- 
space which incorporates information about topology as well 
as morphology of superclusters and voids. ^ 



4 PERCOLATION AND GLOBAL TOPOLOGY 

4.1 Filling factors and one-point function 

A characteristic of the density field which is both simple and 
useful is the one-point probability density function, p(5), 
where 5 = {p — p) /p. In this study we use the integral 



FFc{&) = P{S) = / p{5')d5' 



(7) 



known as the cumulative probability function (cpf) as a ma- 
jor quantity parameterizing the excursion sets. It measures 
the fraction of volume in the excursion set, 5 > 5th. In or- 
der to compare supercluster and void parameters we also 
use the under-density filling factor 



^ Non-geometrical shape-statistics based on mass moments etc. 
can give misleading results when applied to large scale structu re, 
as demonstrated in lSathvaprakash. Sahni fc ShandarinI il99Sl) . 
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Figure 2. Left panel: the fractions of the total volume occupied 
by the largest supercluster (thick dashed line), all superclusters 
but the largest one (thick solid line), largest void (thin dashed 
line), and all voids but the largest one (thin solid line) are shown 
for the density field in the ACDM model smoothed with Ls=5 
/i~^Mpc as a function of the filling factor, FFq. Right panel: the 
y-axis shows the fraction of mass in the components shown in the 
left panel. Vertical dotted lines show the percolation threshold 
for superclusters {FFc ~ 0.07) and voids (FFc ~ 0.78; note: 
FFv = 1 - FFc ~ 0.22). 



Figure 1. Cumulative probability functions of the density con- 
trast in the Virgo simulations of the ACDM model smoothed with 
Ls=5 fc-^and 10 h-^Mpc are shown by the thick and thin solid 
lines respectively. The fractions of mass in the excursion sets are 
shown by the dashed line. The cumulative probability function 
equals the filling factor FFc- 



FFv{5) = 1 - FFc{5) = J p{5')d5'. 



(8) 



The filling factor FFc = FFc{S) is shown in Fig. (thick 
sohd hne) for ACDM smoothed on the length scale of La— 5 
/i~^Mpc. In addition, the thick dashed line shows the frac- 
tion of mass in the excursion set for the same smoothing 
scale. From Fig.Qwe find that the density field is 'nonlinear' 
{5 > 1) in a relatively small fraction ( 15%) of the total 
volume. However, the difference between these two curves 
is the first clear demonstration of nonlinearity of the den- 
sity field. Convergence to the Gaussian distribution with the 
growth of the smoothing scale is demonstrated by the two 
thin lines corresponding to the smoothing scale of Ls — 10 
ft-'Mpc. 



4.2 Percolation 

Understanding percolation is essential for understanding the 
morphology of the supercluster-void network. Percolation is 
important because the properties of superclusters and voids 
radically change at the percolation transitions (see Fig. I^J. 
The left panel shows: (i) the fractional volumes in the largest 
supercluster and void (dashed lines) and (ii) the total vol- 
ume in all superclusters and voids after the largest object has 
been removed from the sample (solid lines) . The right panel 
shows mass fractions in all four components. At relatively 
high thresholds Sc >, 1.8 corresponding to small filling fac- 
tors FFc ^ 0.07 the largest supercluster has insignificant 
volume and mass compared to the total volume or mass 



contained in the over-dense excursion set, 5 > (5th- During 
the percolation transition at FFc ~ 0.07 corresponding to 
5c ~ 1.8, both volume and mass in the largest superclus- 
ter rapidly grow, overtaking the volume and mass in the 
entire excursion set, and completely dominating the entire 
sample from this point onwards. The largest void behaves 
in a qualitatively similar manner if plotted versus FFv- At 
FFv ;< 0.22, Sv < — 0.5 its volume is small compared to 
the volume of the under-dense excursion set, 5 < 5th, but at 
the percolation transition FFv ~ 0.22, Sv ~ —0.5, it takes 
over and from then on remains the dominant structure in 
the under-dense excursion set. Since FFc = 1 — FFv the 
void percolation transition takes place at FFc ~ 0.78 as 
shown in Fig. |21 

Two obvious conclusions can be drawn from the above 
discussion. First, at percolation the object having the largest 
volume becomes very different from all remaining objects, 
therefore it must be studied separately. Second, individual 
objects - both superclusters and voids - must be studied in 
the corresponding phase before percolation occurs in the cor- 
responding phase. Both superclusters and voids reach their 
largest sizes, volumes and masses just before percolation sets 
in. 

Figure |5| also shows that at FFc > 0.5 superclus- 
ters dominate by volume while at FFc < 0.5 voids dom- 
inate. In the range 0.2 FFc ;^ 0.7 corresponding to 
^■^ Si ^ Si ~ 0.43, no more than 10% of volume is occupied 
by non-percolating objects, while the remaining more than 
90% of entire volume is occupied by just two largest objects: 
percolating supercluster and the percolating void. In the 
range between the two percolation thresholds FFc ~ 0.07 
and FFc ~ 0.78, both largest objects percolate and there- 
fore the density field has a sponge like structure. The inter- 
val of sponge like structure in a Gaussian field is between 
FFc ~ 0.16 and FFc ~ 0.84. Therefore nonlinear grav- 
itational evolution has shifted it toward smaller FFc and 
increased its range by a little bit from 0.68 (=0.84—0.16) to 
0.71 (=0.78-0.07). 

It is worth stressing that the shift and length of the 
sponge like interval are determined by the percolation 
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thresholds which in general are two independent parame- 
ters. As a result of these shifts the interval of the so called 
"meat-ball" topology (when no supercluster percolates) has 
reduced compared to the Gaussian case from FFc ~ 0.16 
to FFc ~ 0.07 and the interval of the "bubble" topology 
(when no void percolates) is increased from FFy ~ 0.16 
to FFv ~ 0.22. All the numbers obviously depend on the 
adopted smoothing scale but the sign and type of change 
must be universal for the ACDM model. 

For both superclusters and voids FF — FFmax + J'^ind, 
where Ffind stands for the fractional volume occupied 
by all objects excluding the largest one, and FFm^x is 
the fractional volume in the largest object. Further since, 
i^^max/f'-Find ^ 1 iu the the most part of the range be- 
tween two percolation thresholds for both superclusters and 
voids, it is not surprising that FF^^^ = FF^^^ almost ex- 
actly at FFc ~ FFv ~ 0.5. A similar r esult was found in 
the m odels with power law initial spectra ijYess fc Shandariiil 
\im(t) . 

Percolation is characterized by many features, the most 
conspicuous being the rapid merger of disjoint parts of the 
excursion set into one connected structure spanning the 
entire volume. Merging of superclusters occurs when the 
density threshold is reduced whereas merging of voids cor- 
respond to the growth of the threshold. Spanning of the 
largest supercluster or void throughout the whole volume 
results in connection of the opposite faces of the cubic vol- 
ume by this structure which explains the term percolation. 
Although in principle the percolation transition can be de- 
termined by checking if the opposite faces of the cube are 
connected, in practice it is more robust to identify perco- 
lation using other properties of the excursion set iKlvpinl 
Il987l: iKlvpin fc Shandarinlfl993ll . As the estimators of the 
percolation threshold we use the following four ratios 



MF, 



mf: 



(i) 



1,2,3,4, 



(9) 



where MF^'^ is one of four quantities: volume, area, inte- 
grated mean curvature, or mass 

A, MF^^^ = C, Mf = M). The subscript "max" is 
self-explanatory, and "i?S" stands for the entire excursion 
set. At percolation these ratios grow extremely rapidly from 
very small values to unity. The maximum rate of growth 
can be used as a reliable estimator of the per- 
colation threshold as shown in Fig. |3 All four parameters 
detect the percolation transitions at FFc ~ 0.07 for su- 
perclusters and at FFv ~ 0.22 for voids. These transitions 
correspond to the density contrast Sc ~ 1.8 for superclusters 
and Sv ~ —0.5 for voids. 

Figure 01 compares the largest supercluster with the 
largest void. Three Minkowski functionals and the frac- 
tion of the mass in the largest objects are shown as a 
function of the corresponding filling factor. The difference 
of two curves in every panel is a significant indication of 
non-Gaussianity of the density field. All Gaussian curves 
must coincide in every panel and show maximum growth at 
FFc'"' ~ FPv ~ 0.16. The nonlinear gravitational evolu- 
tion in the ACDM model shifts the percolation transition 
in the over-dense excursion set toward smaller filling factors 
{FFc ■ 0.16 ^ 0.07) and in the under-dense excursion set 
toward larger filling factors {FFv ■ 0.16 ^ 0.22). Again 
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Figure 3. Estimates of the percolation thresholds. The rate of 
growth (5m(*) /5{FF) for the four estimators listed in eq. |5] is 
shown as a function of FF. Thick lines show results for super- 
clusters. All four curves consistently peak at FF = FFc = 0.07. 
Thin lines show similar quantities for voids with a distinct peak 
at FF = FFv = 0.22. Solid, short dashed, long dashed, and 
dot-dashed lines show the volume, area, curvature, and mass es- 
timators respectively. Vertical dotted lines mark the percolation 
thresholds. 
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Figure 4. The mass fraction and Minkowski functionals of the 
largest (by volume) supercluster and void in the ACDM model 
smoothed with the Ls=5 h~^ Mpc window as a function of cor- 
responding filling factor {FF = FFc foi" the largest supercluster; 
FF = FFv for the largest void). Solid and dashed lines show 
the parameters of the largest supercluster and void respectively. 
Vertical dotted lines mark the percolation transitions. 
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the particular numbers depend on the choice of the smooth- 
ing scale but the sign of the effects is independent of the 
smoothing scale. At smaller smoothing scales as well as for 
any adaptive smoothing having a better resolution in high 
density regions the supercluster percolation threshold must 
decrease {FFc < 0.07) and void percolation threshold must 
increase {FFy > 0.22). On the other hand, increasing the 
smoothing scale would result in a continuous reduction of 
differences between superclusters and voids and their ulti- 
mate convergence to Gaussian curves in every panel (not 
shown) . 

ft is worth noting that the percolation transition in 
voids is more conspicuous than that in superclusters. The 
transition is particularly clearly marked by Amax and Cmax 
curves. All curves look differently after percolation as well. 
However, in order to precisely evaluate the significance level 
of these differences one needs to analyze more than one re- 
alization and/or have larger simulation volume. 

4.3 Global topology 

It is interesting to compare the percolation and genus statis- 
tics. Both were suggested as tests for assessing the con- 
nectedness of the large-scale structure. In a series of pa- 
pers Zel'dovich and Shandarin fe cldovich 1982; Shandari3 
Il983l: IShandarin fc ZeldovichI [l983D raised the question of 
topology of large-scale structure and suggested percolation 
statistic as a discriminator between models. The percola- 
tion test was first applied to a rodshift catalog (compiled by 
J. If uchra) by Zcld ovich, Einasto & Shandarin ( f 982^) and 
then lEinasto et al.l lll984lh who found that the connectivity 
between galaxies in this catalog was significantly stronger 
than for a Poisson distribution. In contrast, a non-dynamical 
computer model having approximately correct correlation 
functions up to the fourth order ( Sone ira & Peebles 1978) 
showed significantly weaker connectivity than observed. 
Thus, percolation was able to detect connectedness in the 
galaxy distribution. It was also demonstrated that three low- 
est order correlation functions (two-, three- and four-point 
functions) are not sufficient to detect the connectedness in 
the galaxy distributio n. 

A few years la ter lOott. Melott fc DickinsonI l)l98fitl (see 
also iMelotd l)l99(t) for a review) suggested the genus statis- 
tic as a discriminator between various models of large-scale 
structure. Although both percolation and genus statistics 
characterize density fields and are sensitive to the connect- 
edness of the large scale structure, each carries significantly 
different information. It is important to remember that the 
genus refers to a surface which is the interface between over- 
dense and under-dense regions (defined at a given density 
threshold). An interpretation of genus as a characteristic of 
a three dimensional object (supercluster or void) bounded 
by this surface is not unique and certainly non-trivial if the 
region has a complicated shape. For instance, both a full 
sphere and a doughnut with a bubble in its body have genus 
of 0. 

Figure |S] shows the genus curve in a slightly unusual 
form (see also Sahin, Sathvaorakash & Shandarin (1997)). 
The half of the curve corresponding to high density thresh- 
olds is shown as a function of the over-density filling factor 
FFc (solid line) while the other half corresponding to low 
thresholds is plotted as a function of the under-density fill- 




Figure 5. Left panel: the global genus is shown as a function of 
the fiUing factor for the density field smoothed with Lg = 5 
Mpc. The half of the curve corresponding to high density thresh- 
olds is plotted as a function of FFq (solid line) while the other 
half corresponding to low density thresholds is plotted as a func- 
tion of FFy (dashed fine). For comparison, the dotted line shows 
the Gaussian genus curve having the same amplitude. The ver- 
tical dotted lines mark two percolation thresholds in the ACDM 
[FF = FFc ~ 0.07 and FF = FFy ~ 0.22) and Gaussian field 
{FF = FFc = FFy 0.16). Right panel: the percolation tran- 
sitions in the same density field as indicated; the genus of the 
largest supercluster (solid line) and largest void (dashed line). 
The vertical dotted lines mark the percolation thresholds similar 
to left panel. 



ing factor FFy (dashed line) . This allows to better illustrate 
the deformations of the curve due to nonlinear effects. The 
Gaussian genus curve is symmetric for positive and negative 
thresholds thus both parts of it overlap in Fig. |^ (dotted 
line). The vertical dashed lines mark three thresholds: the 
supercluster percolation threshold at FF — FFc ~ 0.07, 
the void percolation threshold at FF = FFy fa 0.22 
and the both percolation thresholds in a Gaussian field at 
FF = FFc = FFy ^0.16. 

A marked decrease in the amplitude of the genus curve 
compared to the Gaussian curve at small FF is noticeable 
for both over-dense and under-dense excursion sets. (Small 
FF = high density for superclusters and low density for 
voids.) The global genus curve has no significant features at 
either percolation threshold FFc ~ 0.07 or FFy fa 0.22. 
The right panel of Fig. |^ shows the ratios of the genus of 
the largest object Gmax to the global genus of the excursion 
set Ges for both superclusters (solid line) and voids (dashed 
line) . Both percolation curves shown in the right panel would 
overlap in the case of a Gaussian field (not shown) and 
demonstrate the percolation transition at FF ~ 0.f6. The 
other indicators of the percolation transitions (Eq. |^ are in 
excellent agreement with the right panel of Fig. |3 which can 
be seen by comparing Fig. |S] with Figs. |3] and 0] The split- 
ting of the percolation curves shown in Fig. |3 (right panel) 
as well as in Fig. 2] clearly demonstrates non-Gaussianity of 
the density field. 

After percolation takes place the genus of the largest 
object considerably increases (Fig. (SJ which manifests its 
complex shape. The length and therefore filamentarity of 
the largest object radically changes after percolation takes 
place (see Fig. |S| and |7J and cannot be easily interpreted. 
On the other hand, the breadth and thickness, (and con- 
sequently also the planarity), grow much more gradually 
with the growth of the corresponding filling factor and in 
this sense are similar to non-percolating objects. The mean 
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Figure 6. The length, breadth, thickness, mean density, and 
genus of the largest (by volume) supercluster (solid lines) and 
void (dashed line) in the ACDM model smoothed with Lg = 5 
Mpc are shown as a function of of FFq and FFy respectively. 
Vertical dotted lines show the supercluster (FF = FFc ~ 0.07) 
and void {FF = FFv ~ 0.22) percolation thresholds. 
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Figure 7. Planarity and filamentarity of the largest structures 
as a function of the filling factor. The largest supercluster and 
void are shown by solid and dashed lines respectively. Vertical 
dotted lines show the supercluster [FF = FFq 0.07) and void 
{FF = FFv ~ 0.22) percolation thresholds. 



density of the percolating supercluster at the percolation 
transition is psc ~ 4p so that 5sc ~ 3. The corresponding 
value for the percolating void is pv ~ 0.3 and (5v ~ —0.7. 
Both 5qc and (5v are significantly different from the per- 
colation threshold for superclusters (5c ~ l-8p) and voids 
(Sy = —0.5) respectively (see Figure (SJ. 



5 INDIVIDUAL SUPERCLUSTERS AND 
VOIDS 

In this section we carry out comparative statistical analy- 
sis of physical and geometrical parameters of superclusters 
and voids. We stress that the largest supercluster and the 
largest void both are excluded from this analysis. Before per- 




- 0.5 




Figure 8. Left panel: the number density (in (100 h~^ Mpc)~^ 
units) of the most massive superclusters making 90% of all mass 
contained in non-percolating superclusters as a function of the 
filling factor FFq is shown by solid line. The number density 
of the largest voids making 90% of all volume contained in non- 
percolating voids as a function of the filling factor FFy is shown 
by dashed line. Right panel: the number density of the least 
massive superclusters making 10% of all mass contained in non- 
percolating superclusters as a function of the filling factor FFc 
is shown by solid line. The number density of the smallest voids 
making 10% of all volume contained in non-percolating voids as 
a function of the filling factor FFy is shown by dashed line. Two 
vertical dotted lines mark percolation thresholds for superclusters 
{FFc ~ 0.07) and voids(FFv ~ 0.22) 



eolation transition they both are extreme outliers and after 
percolation they have nothing in common with the other 
objects. The total number of individual objects is shown in 
Fig. IHl There are numerous small objects among both su- 
perclusters and voids which dominate by numbers at every 
threshold (right panel of Fig. |SJ|. Including them into the 
statistical analysis along with large objects would seriously 
affect all the parameters. One way to deal with this problem 
would be the computation of weighted parameters, i.e. by 
mass for superclusters and by volume for voids. This ac- 
tually corresponds to computing one or a few moments of 
the probability distribution function which also may be mis- 
leading because the distribution functions are strongly non- 
Gaussian. We analyze only the the most massive superclus- 
ters making 90% of all mass in non-percolating over-dense 
objects and most voluminous voids making 90% of total vol- 
ume in non-percolating under-dense objects. The smallest 
over-dense and under-dense objects are excluded from the 
analysis. Figure |H| shows the number density in (100 h~^ 
Mpc)~^ units for both large (left panel) and small (right 
panel) objects. Despite the smallest objects make only 10% 
of mass or volume they are more numerous then the largest 
objects making the most of mass or volume. 



5.1 Masses, volumes and mean densities 

We begin with the analysis of masses, volumes and mean 
densities 

= ^ 

pmean = "^J" 

which are probably the most important factors among those 
that determine the visual impression. They are also the least 
ambiguous. 

Figure |5] shows the masses, volumes and mean densities 
of superclusters (left) and voids (right) at various thresh- 
olds parameterized by the corresponding filling factor. The 
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Figure 9. Mass (in solar masses), volume (in (h^^ Mpc)^), and 
mean density (i.e. Pmean — (M/V)/p) of all but the largest 
structure are shown as a function of the volume filling factor FFc ■ 
Thick solid line is the median of the distribution, dashed lines 
mark 25-75% interval, and thin solid lines show the third largest 
and third smallest value in the distribution. Vertical dotted lines 
show the percolation thresholds in over-dense and under-dense 
excursion sets. The supercluster and void parameters are shown 
in the right and left panels respectively. 



thick solid lines show the median of the distribution, the 
dashed lines show the 50% interval (25% largest and 25% 
smallest are beyond this interval), and thin solid lines show 
the 95% interval. The intervals are highly asymmetric with 
respect to the median which indicates the non-Gaussianity 
of all distribution functions. Every distribution function has 
a long tail at large values. Both superclusters and voids have 
largest masses and volumes around the corresponding per- 
colation threshold although the maximum is more distinct 
in the case of voids. As one might anticipate, voids are less 
massive and more voluminous than superclusters. However, 
the difference is not huge at the chosen smoothing scale. 
The median mass reaches a maximum about lO^'^ Mq for 
superclusters and lO^'*'* Mq for voids. The median volume 
in the corresponding maximum is about 3 X 10^ {h.-'^Mpcf 
for superclusters and about 3 x 10'' {h~^ Mpcf' for voids. 
The mean densities of the superclusters are well above the 
mean density before percolation and gradually decrease after 
percolation takes place however remaining above the mean 
density in the universe. As expected the mean density of 
voids increases with the growth of FFv corresponding to 
the growth of the density threshold. It is somewhat surpris- 
ing that the volumes of superclusters do not change much 
with the threshold, the masses change more but not greater 
than by an order of magnitude. 
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Figure 10. The length L, breadth B, and thickness T of all but 
the largest object is shown as a function of the corresponding 
filling factor. The notations are as in Fig. |5] 



5.2 Sizes and Shapes 

Three characteristic sizes and shapes of superclusters and 
voids can be estimated from Minkowski functionals of ev- 
ery object (eq. |S] and |5J|. Figure [Till shows the median and 
50% and 95% intervals for the length, breadth and thick- 
ness of the superclusters and voids. The sizes of superclus- 
ters are shown as a function of FFc, while the sizes of voids 
as a function of FFv- It is surprising that the median thick- 
ness of superclusters depends on the threshold so weakly; 
it is within 4—6 h~^Mpc interval for a range of thresholds 
between < 5 < 6. This may indicate that the actual thick- 
ness of superclusters is significantly smaller and the mea- 
sured values reflect the width of the smoothing window. 
The breadth of superclusters is not much larger than the 
thickness and it is likely that this quantity is also affected 
by the width of the filtering window. Voids are a little fatter 
than superclusters and their median thickness reaches about 
9 h~^ Mpc at the percolation threshold. Interestingly voids 
are also wider and longer than superclusters (please note 
the logarithmic scale used for length). Recalling that the 
size parameters are normalized to the radius of the sphere 
rather than to diameter we conclude that the longest 25% 
of superclusters are longer than about 50 h~^ Mpc and 25% 
of voids are longer than about 60 h~^ Mpc. 

The shape parameters of superclusters and voids are 
shown in Fig. 1111 The median planarities of superclusters are 
small which means that in the ACDM universe the dark mat- 
ter density field smoothed with 5h~'^Mpc typically has no 
pancake-like superclusters with the diameters greater than 
about 10 h^^ Mpc. The outliers can reach planarity P ~ 0.3 
corresponding to the ratio B/T ~ 1.4 which is not large ei- 
ther. In addition, it happens only at quite large filling factors 
where the density threshold is quite low and only very few 
superclusters are left (see Fig. |HJ . Voids show an opposite 
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Figure 11. The planarity, P and filamentarity, F of all but the 
largest object is shown as a function of the corresponding filling 
factor. The notations are as in Fig. 151 



trend to superclusters: voids are significantly more planar 
than superclusters and the largest planarities in voids oc- 
cur at small filling factors. Filamentarities are significantly 
higher for both over-dense and under-dense objects: median 
values peak at about 0.35 {L/B ~ 2). The outliers could be 
considerably more elongated F > 0.7 (i.e. L/B > 6). As we 
shall see in the next section, the more massive a supercluster 
or more voluminous a void, the greater is its tendency to be 
filamentary. 



5.3 Correlations between morphological and 
physical parameters 

Both superclusters and voids reach their largest sizes near 
their respective percolation thresholds: FFc ~ 0.07 for su- 
perclusters and FFv ~ 0.22 for voids. Figure m\ shows the 
scatter plots of three characteristic sizes (i, B, and T) ver- 
sus mass for superclusters and versus volume for voids. The 
combined plots are made for FFc ~ 0.06, 0.07 and 0.08 
for superclusters and for FFv = 0.21, 0.22 and 0.23 for 
voids. The solid lines show the radius of a sphere having 
the same volume as a given object {R = (3V/47r)i''^). All 
three sizes show a significant correlation with the mass: the 
greater the mass the greater the thickness, breadth and 
length. The thickness and breadth approximately double 
their values and length grows over an order of magnitude 
when mass increases from about lO^*'^ Mq to 10^^'^ Mq. 
Both the thickness and breadth are considerably smaller 
than the radius J? of a sphere having similar volume for 
large superclusters (M ^ 10^^ Mq) as well as for large voids 
(V > 10" (h~^ Mpcf). On the other hand the length is con- 
siderably greater than R. This is another manifestation of 
anisotropics of the large-scale objects. 

The corresponding plots of planarities and filamentar- 
ities show a similar correlation: the larger the mass of a 
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Figure 12. The length, breadth, and thickness versus mass for 
superclusters and versus volume for voids at percolation. Solid 
circles show the relation at percolation thresholds: FFq = 0.07 
for superclusters and FFy = 0.22 for voids. Crosses show the pa- 
rameters before percolation (FFq = 0.06 for superclusters and 
FFv = 0.21 for voids) and empty triangles after percolation 
{FFc = 0.08 for superclusters and FFv = 0.23 for voids). SoUd 
lines show the radius of the sphere having the same volume as 
the corresponding object. Note the logarithmic scale used for the 
length. Three lines correspond to three different thresholds. 
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Figure 13. The planarity and filamentarity vs mass (for super- 
clusters) and vs volume (for voids) at percolation. Notations are 
as in FiglT^ 
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Figure 14. Right panel: correlation of tlie mean density contrast 
{^mean = M/ (pV) — 1) of a supercluster witli its mass at perco- 
lation. Left panel: correlation of the mean density contrast of a 
void with its volume at percolation. Notations arc as in Fig |12l 
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Figure 15. The maximum of genus as a function of the filling 
factor. Left panel: the maximum genus in isolated superclusters. 
Right panel: the maximum genus in isolated voids. Note the dif- 
ference of the scales on the vertical axis. 



supercluster or the larger the volume of a void, the greater 
its planarity and filamentarity (Fig. . One can clearly see 
that the largest objects (superclusters with M > 10^^ Mq 
and voids with V > 10* {h~^ Mpcf' are the most anisotropic 
large-scale objects. Note that voids display a higher level of 
planarity when compared to superclusters. Indeed, one of 
the most noticeable results of this analysis is evidence that 
voids defined near the onset of percolation of the under- 
dense excursion set are significantly non-spherical. Finally, 
the larger the mass of a supercluster the greater its mean 
density and the more voluminous a void the lower its mean 
density (Fig.[HJ. 



6 SUBSTRUCTURE IN SUPERCLUSTERS 
AND VOIDS 

Visual inspection shows that superclusters and voids 
have noticeable substructure. Isolated voids inside su- 
perclusters and isolated clusters and superclusters inside 
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120031) show that voids are also filled with smaller fila- 
ments of high density. In the ACDM universe these fila- 
ments are strong enough to survive smoothing on a scale 
Ls = 5 h~^Mpc. Figtire lT^ shows the maximum of the genus 
of superclusters and voids. In order to interpret it properly 
we need to explain how to interpret the genus of an object. 
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Figure 16. The genus-mass relation is shown for isolated super- 
clusters in the left hand side panel. The genus-volume relation 
is shown for voids on the right. Every isolated supercluster and 
void having genus greater than zero at all density thresholds is 
shown. The mass is given in the solar units and the volume in 
(h~^ Mpc)^. Note the difference of the scales on the vertical axis. 



A stand ard interpretat ion of genus in cosmological literature 
(see e.Q. lMelotd iimd) ) says 

G — (number of holes)-(number of isolated regions)-f 1, (11) 

this requires additional clarifications. First, "holes" stand 
for fairly complex mathematical objects that do not always 
coincide with the simple minded meaning of the word. For 
instance, a cavity in a tooth is not a hole, a bubble in a glass 
of champagne is not a hole either, but a tunnel through a 
mountain is a hole. Second, "number of isolated regions" 
means the number of isolated pieces of the surface that 
define the boundaries of the object. For instance, a region 
bounded by two nested spheres (a thick spherical shell) has 
two disconnected pieces of the boundary surface and no tun- 
nels, therefore its genus is —1. A torus has surface in one 
piece but it has a tunnel, thus its genus equals +1. Consider 
an object topologically homeomorphic to a sphere (G=0), 
I.e. only one boundary surface and no tunnels. Adding a 
tunnel increases the genus by one and adding a bubble in- 
side the body reduces the genus by one. Thus, a doughnut 
with a bubble of air in its body has genus 0, exactly as a 
full sphere. 

Figure 1151 shows that the genus of an isolated super- 
cluster can be as large as five. It means that there are at 
least five tunnels through the supecluster. The number of 
tunnels could be even greater if it also had a few smaller 
isolated voids inside. The vast majority of superclusters are 
topologically isomorphic to a sphere i.e. , G — 0, but the 
most massive ones have genus greater than unity. A super- 
cluster with genus of unity is homeomorphic to a doughnut 
and one with genus of two to a pretzel. The substructure 
of voids is considerably more complex. The largest genus of 
voids detected in this simulation is 55, therefore at least 55 
filaments span through the void ! 

Figure lTSl shows the correlation between genus and mass 
(for superclusters) or volume (for voids). Largest superclus- 
ters M 10^^ Mq have a nontrivial topology. Generally the 
more massive supercluster the greater chance of complex 
topology. Voids display even stronger correlation between 
genus and volume. 
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7 SYMMETRY OF GLOBAL MINKOWSKI 
FUNCTIONALS 

The globa l Minkowski functionals were int roduced in cos- 
mology bv lMecke. Buchert fc WagneJ lll994 ') and since then 
have been used in many s tudies, see e.g. [Schmalzing ct al. 
Jl999l) : IShethet al.l |20o3). Since the Minkowski function- 
als are additive^ we obtain them by simple summation 
over all objects at every density threshold. As discussed in 
ISheth et alJ ll2003f) ') computing the global Minkowski func- 
tionals poses some problem in the case of periodic boundary 
conditions. However, as shown in Appendix it can be solved 
if Minkowski functionals are computed for every superclus- 
ter and void at a given threshold. 

In Gaussian fields three global Minkowski function- 
als (area A, integrated mean curvature C, and genus G) 
have certain symmetries as functions of the filling factor. In 
addition, the number of positive peaks, A''+ in the excursion 
set equals the number of negative troughs, A'^_ : 



A{1 - FFc) 
C(l ~ FFc) 
G(l - FFc) 
N-{l-FFc) 



A(FFc), 
-C{FFc), 
GiFFc), 
N+(FFc). 



(12) 
(13) 
(14) 
(15) 



The departure from these symmetries manifests the non- 
Gaussianity of the field even before a comparison with the 
corresponding Gaussian functions. Figure [T7I shows the vio- 
lation of the symmetry in the global Minkowski function- 
als in the ACDM model smoothed with Ls=5 h~^Mpc 
and 10 h~^ Mpc. The half of every curve corresponding 
to positive thresholds (iSth >, 0) is plotted as a solid line 
and the other half as a dashed line. The curves corre- 
sponding to smaller smoothing scale have higher amplitudes 
and are shown by thick lines. One can see that while the 
curves corresponding to Ls=5 /i~^Mpc are likely to be non- 
Gaussian,they almost overlap at Ls—lO h~^Mpc. 



8 SUMMARY AND DISCUSSION 

We have studied the large-scale network in the dark matter 
density field in real space obtained from N-bo dy simulations 
by VIRGO consortium i Jenkins et al.|[l993l . The density 
distribution has been smoothed with a Gaussian window 
(L, = 5 h-^ Mpc). 

The major goal of the study was quantitative descrip- 
tion and comparison of over-dense and under-dense regions 
which were dubbed as superclusters and voids respectively. 
We conducted the measurements of four Minkowski func- 
tionals for every individual supercluster and void selected at 
99 density thresholds equispaced in terms of the filling factor 
(Eqs.H|7|and|SJ with the SURFGEN code. By computing 
the shapefinders (Eqs. |S] and [HJ from Minkowski function- 
als we estimated the sizes and shapes of every supercluster 
and void. In this study we did not assume any shape or 
other characteristics for the large-scale objects. Instead we 
have measured the morphology and various other character- 
istics of the objects as they were defined by the isodensity 



^ Note, genus is not additive but the Euler characteristic is ad- 
ditive. 
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Figure 17. Symmetry test for the density field smoothed with 
Ls=5 h-l Mpc (thick lines) and 10 h"^ Mpc (thin lines). Global 
area, curvature and genus at FF < 0.5 are shown as a function 
of 1 — FF by dashed lines. In addition, the sign of the global cur- 
vature is also changed. The number of voids is shown by dashed 
line. The vertical dotted lines mark FF = 0.16 where the Gaus- 
sian genus curve crosses the zero level and the integrated mean 
curvature has a maximum. 



surfaces. It is worth stressing that in our study the super- 
clusters and voids have been defined as three-dimensional 
regions with closed boundary surfaces^. 

The main results of our investigation are summarized 
below: 

• Individual superclusters totally occupy no more than 
about 5% of the total volume and comprise no more than 
20% of mass if the largest {i.e. percolating) supercluster is 
excluded (Fig.l^J. 

• The maximum of the total volume and mass comprised 
by all superclusters except the largest one is reached approx- 
imately at the percolation threshold: 5 « 1.8 corresponding 
to FFc ~ 0.07. 

• Individual voids totally occupy no more than 14% of 
volume and contain no more than 4% of mass if the largest 
void is excluded (Fig.|2l. 

• The maximum of the total volume and mass comprised 
by all voids except the largest one is reached at about 
the void percolation threshold: S « —0.5 corresponding to 
FFv ~ 0.22. 

• Between these two percolation thresholds all superclus- 
ters and voids except the largest ones take up no more than 
about 10% of volume and mass. Both largest supercluster 
and void span throughout the whole space and have a very 
large genus. Therefore they have no well defined sizes, vol- 
umes, masses or easily defined shapes. 

^ In some studies voids are viewed as the regions between fil- 
aments connected in a closed loop; for example, see Fig. 4 in 



ISheth fc Van de Weygaer j <200: 
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• Although superclusters are more massive and voids 
are more voluminous the difference in maximum volumes 
reached at the corresponding percolation thresholds is not 
much greater than an order of magnitude (see Fig. |^ . The 
difference in maximum masses is even smaller than the dif- 
ference in volumes. 

• The volumes, masses and geometrical sizes of super- 
clusters increase as the density threshold decreases and 
reach maximum values at about the percolation threshold 
(5 ~ 1.8). At lower thresholds all parameters decrease as 
the threshold continues to decrease. 

• The sizes of voids are significantly larger than those of 
superclusters even in the density field smoothed with Ls = 
5 h~'^ Mpc. 

• The length of a quarter of the most massive super- 
clusters exceeds 50 h~^ Mpc. The most voluminous voids 
are even longer: 25% of them are longer than 60 h^^ Mpc. 
The longest non-percolating supercluster is as long as 100 
h~^ Mpc and the longest non-percolating void is as long 
as 200 h~^ Mpc. Both are comparable to the size of the 
box (239.5 Mpc) and therefore may be affected by the 
boundaries. 

• The genus value of individual superclusters can be ~ 5 
while the genus of individual voids can reach ~ 55. This 
implies significant amount of substructure in superclusters 
and voids. 

• Voids have considerably more developed substructure 
than supercluste rs. This is in a general agreement with other 
studies of voids (iGrogin fc Gelleill2000t|Peeblesll200j) . 

• One of our main results is that voids, as defined through 
the density field (read dark matter distribution) can be dis- 
tinctly non-spherical. Whether this result carries over to 
voids in galaxy surveys will depend upon the nature of the 
baryon-dark matter biasing and also on whether the density 
field is sampled in real or in redshift space. Since gravita- 
tional lensing probes the density field directly, our results 
are likely to be of some relevance both for ongoing as well 
as future weak lensing surveys of large scale structure. 

• The planarities of both superclusters and voids are quite 
low P < 0.3. This implies that the pancake-like structures in 
the dark matter density in real space are not typical in the 
ACDM model. We are stressing that this conclusion may be 
affected by the size of the smoothing window. 

• The percolation thresholds as well as some other param- 
eters depend on the smoothing scale and for smaller smooth- 
ing scales or adaptive filtering windows the supercluster per- 
colation threshold must decrease [FF^^'^' < 0.07) and the 
void percolation threshold increase {FFy"^'^' > 0.22). 
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10 APPENDIX: GLOBAL MINKOWSKI 
FUNCTIONALS IN A BOX WITH 
PERIODIC BOUNDARY CONDITIONS 

In computing the Minkowski functionals of individual large- 
scale objects (superclusters or voids) we assume that the 
boundary surface of every object is closed however it may 
consist of several disconnected parts. If the region is cut by a 
box boundary it is always closed by the corresponding part 
of the box face. We also assume that the partial Minkowski 
functionals are measured for all superclusters and voids at 
the same set of density thresholds. In this case the evaluation 
of the global Minkowski functionals with periodic boundary 
conditions can be done as follows. 

• Volume: The total volume is obviously independent of 
the assumption of periodicity — l/^*^' = ^ v^'^^ where 
v^'^'' are the volumes of individual superclusters and T/^'^p' 
and V^'^^ denote the global volumes in periodic and non- 
periodic boxes. 

• Area: Let us denote the total area of the superclusters 
at given filling factor FF = FFc as A'^"'' = ^ a*''' and the 



total area of voids as A^"'' = J2 Then, the sum of two 
becomes 



(16) 



where A^^^' is the total area of the excursion set in the box 
with periodic boundary conditions and A^''°^^ is the total 
area of the box boundary (see Fig. 1181 for an illustration) . 
For periodic boundary conditions the global area of the over- 
dense excursion set obviously equals the global area of the 
under-dense excursion set since it is the area of the common 
interface surface. Solving the above equation for one 
obtains 



A 



(17) 



• Integrated mean curvature: The isodensity surfaces is 
constructed by th e triangulation described in detail in 
ISheth et alJ ll2003l) . It is convenient for the further calcu- 
lation of the global integrated mean curvatures to split the 
sum into two parts: one is the sum over the inner edges 
and the other is the sum over the edges lying on the box 
boundary 

all edges in. edges b. edges 

and 

E =^"'= E E (19) 

all edges in. edges b. edges 

For the inner edges the contributions to global mean curva- 
tures of the superclusters and voids have the same magni- 
tude but opposite signs 

1 



(20) 



where ipk is the angle between two surface triangles having 
common edge as shown in Fig. 1181 Please note that in this 
equation we use the angle b etween the triangle planes while 
m eq. 6 jSheth et al.ll2003h we used the angle between the 
normals to the triangles. If the edge lies on the box boundary 
the relation is less trivial. 

Consider a region Ri being cut by the box face into two 
pieces Rn and Ri2 as schematically shown in Fig. 1181 The 
dashed lines show the projections of the opposite faces of the 
box. The plane of the figure is assumed to be orthogonal to 
the edge lying on the faces of the box and having the length 
li. This assumption is made only for convenience of the il- 
lustration and does not affect the calculation. We consider 
the input to the mean integrated curvature from the edge 
li. The code treats the regions Rn and Ri2 as two separate 
regions, therefore 



(21) 



However, if Ri was treated as one region which corresponded 
to the periodic boundary conditions then we would have 



(22) 



Thus, each edge lying on the face of the box contributes to 
the global integrated mean curvature extra ^hn compared 
to the periodic case. Similar reasoning yields for the adjacent 
void 
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Figure 18. Two superclusters are schematically shown. One is 
inside the box (internal) and the other is cut into two pieces by 
the boundary of the box. 



and 



Note that cf^'^ — —c^^' as it should be. 
For the face edges we have 



„(fp) 



(23) 



(24) 



(25) 



Now we have all the terms for equations 1181 and 1191 and af- 
ter simple and straightforward calculations obtain the global 
mean integrated curvatures for superclusters and voids for 
periodic boundary conditions in terms of non-periodic quan- 
tities 

1 



C 



-c 



(vp) 



(26) 



• The Euler characteristic: Three kinds of elements of the 
triangulated surface are triangle faces, edges and vertices. 
The Euler characteristic of the triangulated surface is 



X^F-E + V, 



(27) 



where F, E and V are the total numbers of faces (triangles), 
edges and vertices respectively for both closed and open sur- 
faces. 

A related measure of topology often used in cosmology 
is the genus G = 1 — x/2- Although the two quantities are 
uniquely related and therefore carry exactly the same infor- 
mation the Euler characteristic has more convenient mathe- 
matical properties. As we shall see below including the sur- 
faces with boundaries {i.e. open surfaces) into analysis is 
very useful in many cases in particular in dealing with the 
boundary conditions. For such surfaces the genus is frac- 
tional: e.g. the Euler characteristic of a patch is Xp = 1 and 



thus its genus is Gp — —1/2. As a result the standard in- 
terpretation (eq. IIH cannot be directly applied. Secondly, 
and more importantly, the Euler characteristic is additive 
while genus is not. It means that the Euler characteristic 
of a set of disconnected surfaces is simply the sum over 
all members of the set Xg ~ X/ ^» genus is not 

Gg 7^ '^i9i- example, three isolated spheres have a 
genus of Gg = -2 7^ = while Xa = J2Xi = ^- In 

our study we use Minkowski functionals for both the global 
description of the density field and individual objects. Thus, 
the additivity of a parameter becomes a very useful prop- 
erty. 

Similarly to the previous discussion of the integrated mean 
curvature we split each number into two: one is the input 
from the elements lying inside the box and the other from 
the elements lying on the boundary i.e. faces of the box 



X = {F,n - E,„ + V„,) + {Ft ^Et + Vt), 



(28) 



where subscripts "in" and "b" correspond to the internal 
elements of triangulation lying inside of the box and on the 
boundaries of the box respectively. In the periodic case there 
are neither triangles nor edges or vertices on the faces of the 
box, therefore 



+ Vi„ 



(29) 



for both superclusters and voids because it is the character- 
istic of the common interface surface. 

The boundary group of terms gives different results for 



superclusters and voids. In the case of superclusters F, 



^b^^ + ^b'^' ~ Xb'^' is the Euler characteristic of the "cluster 
boundary surface" i.e. the surface consisting of all regions 
on the faces of the box that have been used for closing the 
superclusters cut by the boundary. And in the case of voids 
it is obviously the remaining of the box surface which can be 
called the "void boundary surface" F^^^ — E^^'' + V^^^ = xi"'- 
For instance, consider a simple case when there is only one 
cluster cut by a box face. Assuming non-periodic boundary 
conditions and closed boundary surfaces for all superclusters 
the parts of cluster surface consists of two simply connected 
pieces each homeomorphous to a circle and thus Xb^^ ~ 2. 
The void boundary surface is homeomorphous to a sphere 
with two holes (not tunnels) and thus x^"' = 0. If there were 



n superclusters cut by the box boundaries then x, 



2n 



and Xj"' = 2 — 2n. In general, since the parts of the boundary 
used for closing the superclusters and the parts used for 
closing voids make the whole surface of the box 



(c) , 

Xb +X 



(v) 



= 2 



(30) 



since the box surface is homeomorphous to a sphere and 
therefore is equal to two. Summing up the cluster ans void 
Euler characteristics one obtains 



X 



+ ^M^2x''>+xr+xl^'^2{x''' + l) 



ip) 



(31) 



Thus, the global Euler characteristic for periodic bound- 
ary conditions, x'^' can be found if both x''^' and x*"' are 
computed at the same threshold under assumption of non- 
periodic boundary conditions 

1 



X 



(p) 



(32) 



Computing x'^'^' and x'"' also assumes that all individual 
clusters and voids have closed surfaces. The surfaces of su- 
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perclusters and voids cut by the boundary of the box are 
closed by necessary parts of the box boundary surface. 
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